Inertial range scaling of scalar flux spectra in uniformly sheared turbulence 
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A model based on two-point closure theory of turbulence is proposed and applied to study the 
Reynolds number dependency of the scalar flux spectra in homogeneous shear flow with a cross- 
stream uniform scalar gradient. For the cross-stream scalar flux, in the inertial range the spectral 
behavior agrees with classical predictions and measurements. The streamwise scalar flux is found 
to be in good agreement with the results of atmospherical measurements. However, both the model 
results and the atmospherical measurements disagree with classical predictions. A detailed analysis 
of the different terms in the evolution equation for the streamwise scalar flux spectrum shows 
that nonlinear contributions are governing the inertial subrange of this spectrum and that these 
contributions are relatively more important than for the cross-stream flux. A new expression for 
the scalar flux spectra is proposed. It allows to unify the description of the components in one 
single expression, leading to a classical K~''^^ inertial range range for the cross-stream component 
and to a new K'^^^"^ scaling for the streamwise component that agrees better with atmospherical 
measurements than the K~'^ prediction of Wyngaard and Cote [Quart. J. R. Met. Soc, 98, 590 
(1972)]. 



I. INTRODUCTION 

Homogeneous turbulence subjected to uniform shear is 
one of the classical problems in turbulence. The under- 
standing and modeling of mixing of a passive scalar in 
such flows is of fundamental interest for the prediction of 
atmospherical flows, dispersion problems and heat trans- 
fer. A key quantity in the mixing and transport process is 
the turbulent scalar flux UiO^ with Ui the velocity fluctu- 
ation and 6 the scalar fluctuation. It is the quantity that 
is accounting for the effects of turbulence on the mean 
scalar profile. In order to characterize the contributions 
of the different turbulent lengthscales to the scalar flux, 
it is important to analyze the spectral distribution of 
scalar flux over wavenumbers, i.e. the scalar flux spec- 
trum. In the literature, the spectral distribution of scalar 
flux produced by a uniform scalar gradient imposed on 
isotropic turbulence, was extensively studied. [1, 2, 3, 4, 5] 
In the present paper attention is focused on the case of 
turbulence in the presence of a constant shear S and a 
cross-stream uniform scalar gradient T. The scalar and 
velocity gradients are chosen in the same direction, that 
we call the z-axis. The mean flow is in the x direction: 



S 



dz 



de 

dz 



(1) 



with Q the mean scalar. The only non-zero components 
of the scalar flux are in this case w9 and ud, the cross- 
stream and streamwise scalar fluxes respectively. The 
three-dimensional spectra corresponding to these quan- 
tities are defined by 



J^u^eiK) = FT/, u,{x)e{x + r) 



(2) 



in which K is the wavevector and FT/j. denotes a Fourier 
transform with respect to the separation vector r. In 



the present paper attention is focused on the spherically 
averaged spectra: 



(3) 



where /j^^^^ ■ (S^{K) stands for integral over a spher- 
ical shell with a radius K, the wavenumber. For large 
Reynolds numbers in the inertial range the cross-stream 
scalar flux spectrum is generally believed to obey: 



dz 



(4) 



as was proposed by Lumley f6l|. e is the spectral energy 
flux. This scaling is identical to the one for the unsheared 
case, in which isotropic turbulence is interacting with a 
scalar gradient (see Bos et al.^ for a discussion of this 
case). 

For the streamwise scalar flux spectrum, Wyngaard 
and Cote 0] proposed the expression: 



Fue{K) ^ ^SK- 
dz 



(5) 



This scaling was however never clearly observed. The at- 
mospheric measurements of Kaimal et al. [8| , Kader and 
Yaglom [oj and Caughey of the streamwise heat flux 
spectrum show an inertial exponent closer to —2.5 than 
to the —3 resulting from the analysis by Wyngaard and 
Cote This discrepancy between the classical scaling 
and observations motivates the present work. 

In section [H] the two-point closure approach described 
in Bos et aL[5j will be extended to the case of uniform 
shear flow. The results of the approach, based on the 
EDQNM (Eddy-Damped Quasi-Normal Markovian) the- 
ory are presented in section Hill and the behavior of the 
spectrum as the Reynolds number is increased is dis- 
cussed. A detailed analysis of the spectra is performed 
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and the scalings observed in the inertial range are shown 
to agree with experimental observations. In section HVl a 
new expression is proposed for the streamwise scalar flux 
spectrum using dimensional arguments. 



II. MODELING THE SCALAR FLUX 
SPECTRUM 

Experimental difficulties limit the study of homoge- 
neous shear flows both in Reynolds number R\ and non 
dimensional time St. Information at higher Rx can be 
obtained from atmospheric measurements in which the 
Reynolds number is typically within the range 10"^ < 
Rx < 10'' (see for example the paper of Bradley et al. 
[ll|). Atmospheric flows, although generally subject to 
shear, can however hardly be considered to be subject 
to a uniform shear and a major problem with atmo- 
spheric experiments is the uncontrolable nature of the 
meteorological conditions. The mean shear is generally 
not constant over a long interval of time, the initial and 
boundary conditions are not easily known by the exper- 
imentalist and the role played by buoyancy effects can 
also be difficult to assess. 

Exactly the opposite is the case with direct numerical 
simulations (DNS): the initial conditions as well as the 
boundary conditions are entirely determined by the scien- 
tist or the engineer. The limits of these 'numerical exper- 
iments' are only dependent on the available computer re- 
sources. These computational limits restricted the DNS 
study of Rogers et al. [l2| to Rx « 40, a Reynolds 
number at which an inertial range is not clearly observ- 
able. Higher resolution DNS were recently performed 
by Brethouwer ^3|, but the results for the scalar flux 
spectra were not presented in this work. The Reynolds 
number limitation can be removed by performing Large 
Eddy Simulations (LES) in which only the large scales 
are resolved: we mention the work of Kaltenbach et al. 
[l^ . The problems are then that the flltering operation 
and the subgrid model can affect the results in the iner- 
tial range. Another problem with both DNS and LES is 
that in homogeneous shear flow the integral lengthscale 
increases monotonically (and much faster than in the case 
of isotropic turbulence) so that, quite rapidly, its size will 
become equal to the one of the computational domain. 

In a recent paper it was shown that the EDQNM 
approach is very suitable to study the Reynolds number 
dependency of the inertial range of the scalar flux spec- 
trum. The study was developed for the case of isotropic 
turbulence. The approach is here extended to homoge- 
neous shear flow. In the case of a fluctuating scalar fleld 
produced by the interaction of an isotropic turbulence 
with a uniform mean scalar gradient, there is an exact 
relation between the 3D spectrum Tuie{K) and its inte- 
gral over a sphere with radius K (Herr et al. [isjV 



:FuAK) {I ~ ^Ji^)FuAK) 



(6) 



ent axis and the wavevector. In the presence of shear, 
relation ^ does not hold anymore; nor the spectral 
tensor can be expressed exactly as a function of 
the wavenumber only. A full EDQNM approach of the 
sheared problem would then require to build and nu- 
merically integrate a wavevector dependent closed set of 
equations [TBI. In order to simplify the numerical task 
that would result of this complete approach, the equa- 
tions are integrated over spherical shells with radius K 
so that the basic quantities in the present approach will 
still be the spherically averaged spectra Fuie{K). The 
difference with the unsheared case will be that the shell 
averaged description is now a simpliflcation that will re- 
quire more closure assumptions in the model derivation. 

The main quantities in this paper are then the spher- 
ically averaged spectra: Fu^eiK) defined by (O and 
ipij{K) defined as: 



JT.(K) 



with 



/t 



{x)uj{x-\-r) 



(7) 



(8) 



The approach consisting in considering and modeling 
only the spherically averaged spectra was introduced by 
Cambon et al. 17] for the velocity field. The averaging 
procedure introduces unclosed terms in the equations of 
both ifij{K) and Fu^eiK). The unclosed terms that are 
related to the interaction with the mean velocity field 
(linear transfer and rapid pressure terms) are modeled 
by tensor invariant theory. The nonlinear terms are mod- 
eled by the EDQNM theory with additional assumptions 
about the anisotropy of the spectra. 

We will calibrate the resulting model by comparison 
with experimental results on the turbulent Prandtl num- 
ber and apply it to flows with Reynolds numbers rang- 
ing from low values, attainable by DNS and laboratory 
experiments, up to values that correspond to the high- 
est Reynolds numbers observed in atmospheric measure- 
ments. We note that the molecular Prandtl number in 
this paper is taken equal to unity. 



A. The equation for the scalar flux spectrum 

In homogeneously sheared turbulence in the presence 
of a uniform mean scalar gradient, the equation for the 
3D scalar flux spectrum is 



dU, 



dUn dKnTu 

dxi dKi 



dxj 



with /i the cosine of the angle between the scalar gradi- 



KiKnK- 



(9) 
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In this equation J^ujO, ^ij and the two-point triple cor- 
relations Tfl,„ , 



Tgin = FT/j.ie{x)ui{x)un{x + r) 



(10) 



are functions of the wavevector and time, and a are 
the molecular viscosity and diffusivity respectively. By 
integration of equation ^ over spherical shells of radius 
K, the equation for is obtained 



de 



+T^,\K,t)+n^^{K,t). (11) 



NLi 



The terms on the left hand side of this equation do not 
require modeling: the dissipation and the production by 
shear are closed terms. On the RHS we find a production 
term involving the spherically averaged spectral tensor 
Lpij [K, t) defined in expression ([Tj) , and two terms involv- 
ing explicitly the mean shear. The first term in which 
the mean shear appears is the linear transfer 



ie 



dUr. 



T.{K) 



'-d^{K). 



(12) 



This term is directly related to the stretching of material 
elements by the mean velocity gradient. It corresponds to 
a pure transfer in wavenumber space, as its integral over 
K is zero. The second term involving the mean shear is 
the rapid pressure term 



T.(K) 



(13) 



It corresponds to the linear part of the pressure-scalar 
fluctuation correlation. The nonlinear transfer and non- 
linear pressure terms are respectively defined as 



T^e\K, t)= [ iK,, (T,* „ - Te,„) dl](if ), (14) 

JT,(K) 

n^^(K,t)= / 



^^^T.,„dS(if). (15) 



The integral of T^^{K,t) over all wavenumbers is zero 
and this term corresponds therefore to a flux in spectral 
space (cascade term). The nonlinear pressure term is 
a destructive term that tends to decorrelate the scalar 
and velocity fluctuations. These two terms could be 
combined as they both arise from the nonlinear triple 
correlation Tgij . It is however informative to separate the 
two contributions, as it allows a more precise analysis of 
the results, as will be shown in section IIII CI All the 
terms on the RHS of pT|) are to be modeled. In the 
following the dependence on K and t will be omitted 
from the equations. 
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FIG. 1: Velocity spectra ip^^iK) and ipuw{K) for 32 < Rx < 
10^ at St = 0.5. 



B. Production by the mean scalar gradient 

In equation (jlip the production by the mean scalar 
gradient involves the spherically averaged spectral ten- 
sor tpij . Choosing the orientation as given in expression 
([1]) the only relevant components are fww and ipuw ■ The 
spherically averaged spectral tensor is modeled by an ex- 
tension of the EDQNM theory for anisotropic turbulence. 
The extension of this type of closure to anisotropic turbu- 
lence using spherically av erag ed quantities was first pro- 
posed by Cambon et al. In the present approach, 
as in Cambon et al. the equation for ipij is solved. The 
unknown terms appearing in this equation are modeled 
as proposed by Cambon et al. for the pressure contribu- 
tion, and using simpler expressions for the transfer terms 
(Touil et al. [19.]). A detailed discussion of the results 
of this model is beyond the scope of this paper. In the 
present work we only show the relevant spectra {(pww and 
Vuw) for varying Reynolds number. These are presented 
in figure [T] for St — 0.5. It can be observed that, as the 
Reynolds number increases, the slopes of the spectra in 
the inertial range are approaching the classical values of 
—5/3 and —7/3 for the diagonal and non diagonal com- 
ponents respectively. 
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C. Rapid pressure 

The rapid pressure term appearing in the scalar flux 
spectrum equation (eq. Ilip is: 



dxi 



T.(K) 



(16) 



At the level of description of the present model, the inte- 
gral of J-ujO is known, but not the integral of its moment: 



(17) 



In order to express , the tensor Hl^ is modeled assum- 
ing that it can be represented as an isotropic tensorial 
function of Fu^e and ai — Ki/K. This does not imply 
that Hf^, or the resulting model for the rapid pressure, 
are isotropic, but only that their anisotropic character re- 
flects the anisotropy of the F^.g vector. Using the theory 
of invariants and the representation by tensorial isotropic 
functions (as in Eringen [20] or Schiestel |2l|) and retain- 
ing only the lowest order terms, it is found 

+CaianFu^g + D^UiajFu^e + F)'' ajUnFuie (18) 

This expression needs to satisfy symmetry of the indices 
j,rt, incompressibility and the deflnition of F^jO^ expres- 
sion Q. These constraints allow to reduce the number 
of constants from 6 to 1. The resulting expression for the 
rapid pressure term is 



dU. 



n- = {12A- 4)^F„^, - (8A- - 3)^i^„^.,. (19) 



dxi 



D. Linear transfer 

The problem of modeling the linear transfer is strongly 
connected to the modeling of the rapid pressure term. As 
a matter of fact, after some algebra, it is straightforward 
to show that it only requires a modeled expression for 
Hf^ , a problem that has been solved in the previous sec- 
tion. 

The linear transfer is: 



-LiO — 



dxi 



S(A') 



(20) 



This expression can be rewritten as: 



dUn d 



K 



dxj dK Jo ys(A') dKj 



-dY.{K)dK. (21) 



Noticing that the double integral can be rewritten as an 
integral over the volume of a sphere with radius K and 
applying the Gauss divergence theorem to this sphere. 



expression ([2T|) can be rewritten. The result has the form 
of a surface integral involving the outer vector Kj/K: 



T 



L_dUn d 



dxi dK 



K 



(22) 



with defined by (fT7|) . Using the modeled expression 
(fT5)) . the closed form for the linear transfer reads: 



rpl 
J it 



dxj 



dUA d 



dxj I dK 



KFuAK). (23) 



It has to be noticed that the approach used in this paper 
to model the linear transfer is similar to the one proposed 
by Cambon et al. [13, ll^] (see also Clark and Zemach 
|2a |) for the linear transfer of kinetic energy. It can be 
argued that even if expression ([T5)) may not necessarily be 
a good model for the tensor itself, it could still lead to 
a satisfactory representation of the spherically integrated 
expressions ([T6l) and ([20| . It is also stressed that the 
exact forms of the rapid pressure and linear tranfer terms 
do not influence the asymptotic scalings of the spectra in 
the inertial range that are essentially determined by the 
nonlinear interactions. 



E. Nonlinear transfer and nonlinear pressure 

In the presence of shear, the nonlinear transfer and 
pressure term can not be expressed exactly as a function 
of the wavenumber only. A full EDQNM approach of the 
problem would then require to build and numerically in- 
tegrate a wavevector dependent closed set of equations. 
This approach is complex and numerically expensive. We 
will here treat the nonlinear transfer and slow pressure 
terms with the EDQNM model, derived for unsheared 
turbulence, using the formulation of Bos et aZ. Q Apply- 
ing this model, originally developed by Herr et aL[l^ for 
isotropic turbulence, to the case of an anisotropic veloc- 
ity field has to be considered as an approximation. By 
introducing this approximation, a detailed description of 
anisotropy in wavevector-space is lost. The approach is 
justified in the limit of a small shear and is here used for 
strong shear as an extrapolation. The expressions for 
the nonlinear terms T^^ and 11^^ are not reproduced 
here (Equations (14) and (15) of reference 0). These 
closed terms are exactly the same as in reference intro- 
ducing two model constants A' and A". 



III. RESULTS 

A. Calibration 

The model for the scalar flux spectrum contains three 
constants. Two of these constants are the ones involved 
in the EDQNM closure as mentioned in the previous sec- 
tion. These constants were determined in reference 
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and are not changed in the present work (A' = 0.52 
and A" = 0), to satisfy compatibility with reference!^ 
in the case of vanishing shear. The third constant is 
, the one that appears in the hnear transfer and rapid 
pressure (expression (fT9|) and ((23)) ). This last constant 
is determined by comparison of the value of the turbu- 
lent Prandtl number with the wind tunnel experiments of 
Tavoularis and Corrsin |24] . In this experiment the mean 
shear is generated by using an array of parallel channels 
of different mean flow with cylindrical heating rods at 
their exits. The shear rate in the simulations is cho- 
sen equal to the one in the experiment, S = 47s~^. The 
turbulent velocity field, assumed to be initially isotropic, 
is initialized by a von Kdrmdn spectrum. [25| This ini- 
tial spectrum is chosen so that the Reynolds number at 
St = 12 is R\ ~ 150, corresponding to the experimen- 
tal values (130 < R\ < 160). The scalar fluctuation 
is initially zero. Nonzero scalar fluctuations would not 
influence the scalar fluxes as long as they are initially 
isotropically distributed. 

The turbulent Prandtl number is the ratio of the 
Reynolds stress to the cross-stream turbulent scalar flux, 
normalized by the mean gradients: 



r uw 



(24) 



Best agreement for this quantity with the Tavoularis and 
Corrsin [2^ experiment (Pry = 1-1 at St = 12) is ob- 
tained for a value of = 0.385. This value is retained 
in the following. We note that the spectral exponents are 
not sensitive to the choice of this constant. 
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2; Cross-stream and streamwise scalar flux spectra at 
0.5 for 32 < 7?A < 10*. 



B. Spectra and Reynolds number dependency 

In flgure [2l the spectra of the streamwise and cross- 
stream scalar flux are plotted at St — 0.5, for six 
Reynolds numbers in the range 32 < i?^ < 10^. The 
streamwise component Fug{K), shows a steeper inertial 
range than F^g [K] . Figure [3] shows the spectral expo- 
nents of these two components of the scalar flux as a func- 
tion of the Reynolds number. The exponents are deter- 
mined at two different times, corresponding to St = 0.5 
and 5^ = 12. It is observed that the difference between 
the slopes at different times decreases with increasing 
Reynolds number. The spectral exponent tends faster 
to its asymptotic value for St = 0.5 than for St — 12. 
A possible explanation is the following: the nonlinear 
triple correlations that are responsible for the creation 
of an inertial range are initially zero. This yields an ini- 
tially weak nonlinear transfer, corresponding to a steep 
spectrum. When the triple correlations are fully devel- 
oped, the nonlinear transfer becomes stronger, resulting 
in a less steep spectral distribution. The cross-stream 
scalar flux spectrum tends to the Uyjg = — 7/3 asymptote 
for large R\, like in an isotropic turbulence with a mean 
scalar gradient. The asymptotic value is not changed 
by the presence of shear, as was already anticipated by 



Lumley 6]. As is the case in the absence of shear, the 
asymptote is approached very slowly when increasing the 
Reynolds number (see the recent paper by the authors 

:5]). 

The streamwise spectrum (fig. [51 bottom) is found 
to be steeper than the cross-stream spectrum. Its expo- 
nent, n„e (fig. [3l bottom) tends to a value significantly 
larger than 7/3. However, the high Reynolds number 
asymptote is found to be smaller than the value Uue = 3 
proposed by Wyngaard and Cote 0, and n^g = 23/9 ap- 
pears a more plausible value. This can also be observed 
in figure IH where the spectra compensated by if^^/^ are 
plotted. The normalization used in this figure will be 
clarified in section pvp . 

It has to be noticed that the spectra shown and ana- 
lyzed in this section were obtained at moderate or large 
nondimensional times, St = 0.5 and St = 12. To analyze 
the behavior at shorter time, the time evolution of the 
streamwise scalar flux spectrum is plotted in figure [3 It 
can be observed that at very short times, the spectrum 
behaves as K^^l^ ^ a scaling in agreement with the obser- 
vations of Antonia and Zhu [2^ , who showed two cases of 
spectra with a K~^l^ scaling for over two decades. The 
present model therefore can produce a scaling in agree- 
ment with the one proposed by these authors. However 
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FIG. 3; Inertial range slopes of the scalar flux spectra as a 
function of the Reynolds number. Top: cross-stream scalar 
flux, the dotted line corresponds to K~''^^ . Bottom: stream- 
wise scalar flux, the dotted lines correspond to K~^^^^ and 



this scaling is only observed at very short time and does 
not persist for a longer period and rapidly the spectra be- 
come steeper, tending to the asymptotical if~23/9 bej^ay- 
ior. We note that, starting from an initial energy spec- 
trum with a K~^^^ distribution, the result n^g — —5/3 
can be obtained by purely linear arguments, neglecting 
all the nonlinear terms. 



C. Budget of the different terms in the scalar flux 
spectrum equation 

The different contributions to the evolution equation 
of the scalar flux spectra, equation (fTTj) . are shown in fig- 
urc[6]at R\ ~ 2384. To facilitate the analysis and reduce 
the number of terms, the contributions of the produc- 
tion by the scalar gradient and mean shear were lumped 
together: 



Pi = -Tipis - S6iiF^ 



(25) 

It can be directly deduced from that there is no 
explicit contribution from the mean shear to the pro- 
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FIG. 4: Compensated cross-stream (top) and streamwise 
(bottom) scalar flux spectra for 32 < i?A < 10* at St = 0.5. 



duction term of F^e- As for the dissipative term, 
— {v + a)K'^Fuie, it is negligibly small at this R\ and 
is therefore not shown. 

For simplicity the indices are dropped (all contribu- 
tions to i^„e should have an index 1, the contributions to 
F^g an index 3). All the contributions to F^g, which is a 
negative spectrum, are multiplied by — 1 to facilitate the 
comparison with F^g , which is a positive spectrum. 

It can be noted that for both spectra F^g and F^g , the 
main contributions are the production P and the nonlin- 
ear transfer and pressure terms, and 11^^ respec- 
tively. The linear transfer and rapid pressure II'^ 
terms are relatively small, especially at high wavenum- 
bers. By comparison with Bos et aL[5|, in which the case 
of isotropic turbulence with a uniform scalar gradient was 
analyzed, it is observed that the behavior of P, 11^^ and 
rpNL roughly the same. It can be concluded that the 
presence of the mean shear does not drastically affect the 
mechanisms involved in the dynamics of F^g. Compared 
to the isotropic case, the presence of the mean shear does 
not directly affect F^g, as its production remains by the 
mean scalar gradient only. F^g is essentially produced 
at the large scales and destroyed by the pressure effect. 
This destruction takes place locally as well as at smaller 
scales, since the transfer term introduces a cascade pro- 
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FIG. 5: The streamwise scalar flux spectrum at four diflferent 
times in the beginning of the calculation. The spectra are 
compensated by K^^^ . 



cess. 

The streamwise scalar flux spectrum F^e does not exist 
in the isotropic case. It can be observed in figure [6] (bot- 
tom) that Fuo is essentially generated at large scales by 
the production term associated with the presence of the 
mean gradients. This was already the case for Fu]e, but 
for Fue both the mean shear S and the scalar gradient 
r contribute to the source term. The transfer term still 
introduces a cascade effect and the destruction remains 
associated with the nonlinear pressure. The main differ- 
ence with the Fw0 budget is that in the case of F^e , the 
mean production falls off more rapidly as the wavenum- 
ber increases. The two contributions to Pi in equation 
(pSj) indeed behave as K~'^/^, whereas in the F^e case 
they were found to scale as K~^/^. As a result of this 
more rapid extinction of the production term, it can be 
observed in figure[6l that in the inertial range of the spec- 
tra, there are only two dominant terms: the nonlinear 
transfer T^^ and the nonlinear pressure term 11^^. 



IV. DIMENSIONAL ANALYSIS 

In this section, a dimensional analysis compatible with 
the results of the spectral closure is proposed. First, the 
classical dimensional analysis is recalled, and it is stressed 
that it leads to a spectral scaling that is not compatible 
with the results of the previous section. Then, based on 
the observations in section III.c, a new argument is pro- 
posed. It leads to a scaling that agrees with the EDQNM 
results, as well as with existing atmospherical measure- 
ments. 

For the streamwise scalar flux spectrum dimensional 
analysis based on the quantities S, e and K provides the 
following expression for the spectrum: 
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FIG. 6: Contributions to the equations of the cross-stream 
and streamwise scalar flux spectra at St = 12 and Rx = 2384. 



This expression is linear in the scalar gradient, as it 
has to be to reflect the linearity of the scalar equation. 
Linearity in S is not mandatory since the Navier-Stokes 
equations are not linear; if linearity in S is assumed, (|26|) 
reduces to the Wyngaard and Cote formulation: 



BP) 

FueiK) ^ ^SK-\ 
az 



(27) 



This formulation can also be found by assuming an equi- 
librium between the production term scaling as 



oz 



(28) 



and the pressure destruction term assumed to scale as: 

1 



-Fue{K), 



(29) 



S°'e—K- 



(26) 



tnl{KY 

in which tnl{K) is the nonlinear cascade time 

rNL{K) ~ 6-V3ir-2/3. (30) 

The analysis of the different contributions in figurelHlhav- 
ing shown that the nonlinear terms are dominant in the 
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inertial range and that the local production term plays 
a negligible role, an alternative reasonin g, b ased on a 
wavenumber dependent flux (cf. Bos et a/. [23]) could be 
defendable. For this reason it is proposed to introduce a 
tensorial extension of ([4]): 



-7/3 



(31) 



in which (K) is a quantity that reflects the anisotropic 
character of the nonlinear cascade and has the dimension 
of e^^^. We define A such that each of its components 
is the cube root of the corresponding component of : 



A,: 



1/3 



in which stands for a tensorial spectral 



flux, which is not necessarily conserved in the cascade: 
eij = Cij (K) . It can be assumed that Cij (K) is a spectral 
anisotropic flux in wavenumber space, whose anisotropy 
will be proportional to the anisotropy of the spectral ten- 
sor iptj{K): 



E{K) ' 



(32) 



in which e represents the total energy flux and E{K) 
is the turbulent kinetic energy spectrum. Using Kol- 
mogorov scaling for the energy spectrum and Lumley's 
expression 



(33) 



for the Reynolds-stress spectrum, expression yields 
in the inertial range: 



tww{K) = e. 



(34) 



which yields the Lumley scaling (jl]) for the cross-stream 
scalar flux and for Fue{K): 



CONCLUSION 



Two point closure theory was applied to study the 
scalar flux in homogeneously sheared turbulence with a 
cross-stream scalar gradient. The equation for the spec- 
trum FuioiK) was closed by EDQNM theory and tensor 
invariant modeling. The resulting model was used to 
study the Reynolds number dependency of the spectra 
associated to the cross-stream and streamwise scalar flux. 
The asymptotic value of the spectral exponent of the 
cross-stream scalar flux spectrum is shown to be —7/3, 
similar to the case without shear. The spectral exponent 
of the streamwise scalar flux spectrum is shown to ap- 
proach a value close to a —23/9 scaling, disagreeing with 
classical dimensional analysis which predicts a —3 scal- 
ing. The —23/9 scaling appears to be in better agreement 
with atmospheric observations. 

A new dimensional argument was then proposed, de- 
pending on an anisotropic spectral energy flux. It leads 
to an expression in agreement with both the ii'~23/9 |^g_ 
havior of the cross-stream component and a K~'^^^ de- 
pendency for the streamwise spectrum. 

It is sometimes considered that two-point closure the- 
ories can only yield trivial results, entirely predictable 
by dimensional arguments. The possible existence of the 
^-23/9 scaling illustrates the opposite. Indeed, this scal- 
ing can be obtained by dimensional analysis as shown in 
section IIVI but a K^'^ scaling can also be obtained by 
dimensional analysis. The authors hope that the present 
two-point closure results might inspire future work, which 
prove or disprove the scaling 



F^o{K) ^ ^5i/3e2/9A'-23/9 (35) 
oz 

This scaling is in full agreement with the results pre- 
sented in this paper. Noticing that 23/9 « 2.55, it is also 
in agreement with atmospheric measurements leading to 
Uue ~ 2.5. We insist here however that the presence of 
the parameter a in (|26[) allows for a wide range of possi- 
ble scalings, so that further research is needed to confirm 
or disqualify the scaling ([35]) . An interesting feature of 
expression ([55)1 is that it is not linear dependent on the 
mean velocity gradient. The implications of this for one- 
point models of the scalar flux deserve more attention. 
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